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A STUDY OF SUPERSONIC SURFACE SOURCES— 

THE FFOWCS WILLIAMS-HAWKINGS EQUATION AND THE 

KIRCHHOFF FORMULA# 


F. Farassat* and Kenneth S. Bren trier' 

NASA Langley Research Center, 

Hampton, Virginia 

Abstract 

In this paper we address the mathematical problem 
of noise generation from high speed moving surfaces. 
The problem we are solving is the linear wave equation 
with sources on a moving surface. The Ffowcs Will- 
iams-Hawkings (FW-H) equation as well as the govern- 
ing equation for deriving the Kirchhoff formula for 
moving surfaces are both this type of partial differential 
equation. We give a new exact solution of this problem 
here in closed form which is valid for subsonic and 
supersonic motion of the surface but it is particularly 
suitable for supersonically moving surfaces. This new 
solution is the simplest of all high speed formulations of 
Langley and is denoted formulation 4 following the tra- 
dition of numbering of our major results for the predic- 
tion of the noise of rotating blades. We show that for a 
smooth surface moving at supersonic speed, our solu- 
tion has only removable singularities.Thus it can be 
used for numerical work. 

1. Introduction 

The problem of noise generation from moving bod- 
ies is very important in aeroacoustics. Two current 
methods of attacking this problem are the acoustic anal- 
ogy and the Kirchhoff formula for moving surfaces. The 
acoustic analogy method is based on different forms of 
the solution of the Ffowcs Williams-Hawkings (FW-H) 
equation 1 . This is a linear wave equation with sources 
on a moving surface. The Kirchhoff formula is also 
derived from a linear wave equation with sources on a 
moving surface 2 . For simplicity, we refer to this wave 
equation here as the Kirchhoff (K) equation. Using 
generalized function theory, both the FW-H and the 
K equations can be written with inhomogeneous source 
terms involving the Dirac delta function with support on 
the moving surface / = 0 and the first derivatives of this 
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delta function. The method of solution of these two 
equations are, thus, identical. 

Obtaining various forms of the solution of these 
equations for subsonic surfaces is fairly easy. These 
solutions (Formulations 1 and 1 A of Farassat ) have been 
published elsewhere 2 " 5 . We will not, therefore, address 
the subsonic case here. We mention that the common 
forms of the solution for subsonic surfaces have a Dop- 
pler singularity which make them unsuitable for super- 
sonically moving surfaces. To obtain new forms of 
solution of the FW-H and the K equations for supersonic 
surfaces, we must integrate the Green’s function of the 
wave equation in a different way than the subsonic 
case 1 ' 2 . This was fully recognized by Ffowcs Williams 
and Hawkings 1 and they laid the foundation for the 
work we present here. The solution of the supersonic 
problem is considerably more difficult than the subsonic 
case and it has taken a lot longer to fully overcome the 
many mathematical obstacles. 

To understand the nature of the complexities 
involved, one must recognize that the problem as treated 
here is four dimensional. We are interested in formula- 
tions which are suitable for efficient numerical noise 
prediction from rotating machinery. This requirement 
puts a restriction on what forms of the solution of the 
FW-H and the K equation are acceptable to us. In 
practice, it has been found that the common formula- 
tions for subsonic surfaces are much more efficient than 
the supersonic formulations even if the latter can also be 
used for subsonic surfaces. Thus, one is forced to use 
more than one formulation in any noise prediction code 
based on the FW-H and K equations. 

In noise calculation, a moving surface, such as a 
blade, is divided into panels and the noise generated by 
each panel is summed up to get the total noise from the 
surface. This means that the FW-H and K equation must 
be solved for an open surface, e.g., a panel on the 
moving surface. Thus we must solve these two equa- 
tions with inhomogeneous source terms that have a 
Heaviside function multiplying the Dirac delta functions 
which describe the open surface. The mathematical 


1 

American Institute of Aeronautics and Astronautics 



treatment of these inhomogeneous source terms pro- 
duces an additional complexity in obtaining solutions of 
these equations for the supersonic case. 

In the subsonic case, the Doppler factor appears as 
the result of integration of the source time variable in the 
Green’s function solution of the wave equation after a 
Lagrangian frame is introduced in which the surface is 
time independent. This step makes the problem essen- 
tially three dimensional. To get other forms of the solu- 
tion suitable for a supersonically moving surface, one 
must use integration over the influence surface of the 
observer space-time variables (jc, t) which we call the 
I - surface . This surface is more fundamental to the 
solution of the wave equation than the actual moving 
surface over which the subsonic formulations are 
integrated. 

In the next section the governing equations of the 
problem under consideration are presented. In Section 3, 
we will give a new solution of the FW-H and the K 
equations in closed form for supersonic surfaces. In Sec- 
tion 4 we will show that the singularities of the solution 
are integrable. The concluding remarks follow. 

2. The Governing Equations 


We next find all the surface contributions of the last 
term of Eq. (1) by taking the space derivatives explicitly 
and using the rules of generalized differentiation 2 * ' 6,7 . 
We get 


□V = ^{[p«„-(p-p 0 )v n ]8(/)} 


d_ 
dx 

dT iJ 


+ aF' [pv » M *- (p - p ° )c n iW) 


+ ^r,8(/) + 


5% 

1 J 


= Q\ + Q 2 + Ql + 04 


(2) 


We will then consider an open surface described by 
/= 0, / > 0 where / = / = 0 is the equation of the 
edge of this open surface 2,6 . We define / such that 
V/ = v where^v is the unit inward geodesic normal to 
the edge / = / = 0 2,6 . To calculate the noise from 
this open surface, we must multiply 8 (f) in Q h Q 2 and 
Q 2 by the Heaviside function H(J ) . We will next use 
the concept of restriction of a variable to the surface/ 35 
0 and then take the derivatives of Q\ and Q 2 terms 
explicitly 6 . We use a tilde under a symbol to signify 
restriction. 


The Ffowcs Williams-Hawkings Equation and the 
Kirchhoff equation are quite well-known in 
aeroacoustics. We will need a special form of these 
equations suitable for our work. The FW-H equation for 
a moving surface f(x, t) = 0 where /> 0 outside the 
body is 


Introducing the notations 


E = pu n -(p-p 0 )v n 

(3a) 

E i = P v « M ,-(P-Po) c2 «t> 

(3b) 


we have 


□V = ^{[pM„-(p-p 0 )v„]8(/)} 

- 4“ {[p(u » - v »)“/ + P n iW)} 



dX;dX; 

1 J 


[ TijH(f)] 


( 1 ) 


2 

where p is (p - p 0 )c , p is the density, and p 0 and c 
are the density and speed of sound of undisturbed 
medium. The local normal fluid and body velocities are 
denoted by u n and v n , respectively. The Lighthill stress 
tensor is denoted and p is the surface pressure on 
/= 0. Note that we assume that the surface /is defined 
such that V/ = n where n is the unit outward normal 
to this surface. The Heaviside function is denoted H(f). 
As proposed by Ffowcs Williams and Hawkings 1 , the 
moving surface /= 0 can be penetrable and we assume 
so here. 


Qi = §- ( lE //(/)§(/)] 

= EH(f)8(f)-E v v 8(/)8(/) 
-Ey n H(f)V(f) (4) 

Q 2 = ^\E. l //(/)§(/)] 

= V 2 • e t +E i v i Hf)Hf) 

-2 H f E n n i H(J)Kf) 

+ E i n j H(fW(f ) ( 5 ) 
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Here, v n is the local velocity of the edge along the geo- 
desic normal v with components v ; , Ej is the projection 
vector E, on the surface of/= 0 and H f is the local mean 
curvature of the surface/= 0 8,9 . The surface divergence 
of Ej is V 2 • E T 9 . Using Eqs. (4) and (5), we get 

01 + 0 2 + 03 = ?i H(/)b(f) + q 2 H(f)b'(f) 


+ q 3 §(/)§(/) (6) 

where we have defined the following symbols 

?1 = d -^>h + E + V 2 E T -2H f E n ( 7 ) 

q 2 = c 2 (p-p 0 )(M^-l) (8) 

?3 = P( v „«v - M « v v) + (P - Po> v » v v ( 9 ) 


1 0OO 8(/) dy = \QdS (15) 

/= o 

where //)is the mean curvature of /= 0 8 ' 9 . Introduce a 
new generalized function (distribution) 8 '(/) by the 
following relation: 

J0(j)§;(/)^ = - J ^ds (16) 

/= 0 

The subscript 5 in 8 '(/) stands for “simple” which 
emphasizes the similarity of 8'(/) to the one dimen- 
sional S'(x) that behaves as follows: 

J <|)(x) S'(x) dx = -f(0) (17) 

Now using the results of eqs. (15) and (16) in Eq. (14), 
we see that the following relation holds: 


In Eq.Eq.Eq. (7), E„ = £)«,■ and M n = vjc is the local 
Mach number on / = 0. Note that q 2 in Eq. (6) is 
restricted to the surface /= 0. Also note that in Eq. (5), 
we have dropped the restriction on any variable that 
multiplies 8 (f) if it is not differentiatedPhus we write E 
and not E . There is also no need to use restriction sign 
on V 2 • E t since E T is already restricted to /= 0. It is 
important to recognize that E is the rate of change of E 
as measured by an observer on the surface. 

The FW-H equation for an open penetrable surface 
moving at supersonic speed is: 

□V = ?i H(fMf) + q 2 H{f)h\f) 

+ q 3 H/)Kf) (10) 


8'(/) = 8;(/) + 2/7 / 8(/) (18) 

Equation (18) is next used in Eq. (10) which is writ- 
ten as 

□y = q x H(f) 8 (/) + q 2 H(J) 8 '(/) 

+ ^3 8 (/) 5 (/) (19) 

where now only the definition of q^ is changed as fol- 
lows: 

c)T u 

q \ = n i + E + SJ 2 -E T 

+ 2H f [pv n (v n - M„)-p 0 v„ 2 ] (FW-Heq.) (20) 


A similar equation is also obtained for derivation of the 
Kirchhoff formula 2 for an open surface/= 0, / > 0 

where we have: 


<h 


y_i M ^ 

dn c " dt 


l -^(M n p') + 2H f p' (11) 


q 2 = (M 2 n -\)p' (12) 

q 3 = M n M v p' (13) 


and 


dp' 1 *,r dp' 1 d 


+ 2H f M n 


(K eq.) 


( 21 ) 


In the following section, we give the full solution of 
Eq. (19). 

Remark . We will not address the solution of the 
FW-H equation with the pure quadrupole term alone: 


We will go one further step here in preparation of 
obtaining the new formulation. We note that for a sur- 
face /= 0, IY/1 = 1 , we have the following results 2,6 : 


J 0(.v) §'(/) dy = J 

/= o 


~dQ + 2H f Q 
dn J _ 


dS (14) 


□y-0 4 ^.^(/) (22) 

1 J 

The solution using the collapsing sphere approach is 
singularity free and is given elsewhere 1,6 . 
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3. Solution of Wave Equation With Sources on A 
Moving Open Surface 

In this section, we give the solution to the following 
three wave equations: 


□ 2 ^, = qi (x,t)H(f)d(f ) 

(23) 

□ 2 ^) 2 = q 2 (x,t)H(fW s (f ) 

(24) 

□ 2 4* 3 = ? 3 (*, t) 8(/)8(/) 

(25) 


The source terms here are similar to those of Eq. (19). 
The treatment of these equations are discussed in two 
references by Farassat 2,6 . We will use the solutions to 
eqs. (23) and (25) given in these references but will give 
here a new and particularly simple solution of Eq. (24). 
The materials presented in these references are essential 
in understanding what follows. 


Let F(y; x, t) = f(y, t-r/c) = [f(y,x)] ret , 
and F(y;x,t ) = [f(y, x)\ et , where the sub- 
script ret stands for retarded time. The influence surface 
of the open surface / = 0, / > 0 is called the 
X - surface and is described by F = 0, F > 0 2,6 . The 
edge of this surface is the L-curve described by 
F = F = 0 . Below, we use (x, t) and (y, x) as the 
observer and the source space-time variables, respec- 
tively. The solution of Eq. (23) is 2,6 : 


4ro|> 1 (x, t) = J jA qx \ ret dL 

F = 0 
F> 0 

where 


(26) 


A 2 = 1 +M 2 n -2 M n cos9 (27) 

Here, M„ is the local normal Mach number of the sur- 
face/^ 0 and cos 8 = n r where n is the unit outward 
normal to /= 0, r = (. x-y)/r is the unit radiation vector 
from the source to the observer and r = \x-y\. 

We now consider Eq. (24). The formal solution of 
this equation using the Green’s function method is 

4n$ 2 (x,t) = j^q 2 (y,x) H(f)d' s (f)d(g) dydx (28) 


where g = x - t + r/c. We now introduce a new local 
frame (u\ n 2 , it) where i? = f and u l and u 2 are the 
Gaussian coordinates on /= constant, extended from 
/= 0 along local normal. We assume that w 1 is the length 
variable along the projection of r on the local tangent 
plane and u 2 is the length variable along r x h . Let g (2 ) 


be the determinant of the coefficients of the first funda- 
mental form in the new variables. Let u — > g , then Eq. 
(28) becomes 


r ? 2 CV, TO 3 

4n<\ > 2 (x, 0 = J — - — H( f) 8'(m ) 


x § (g) Jg ( j^du d\? du dx 


-<T 


TO 


H2) 

r sin 8 


Jg = 0 


3 2 3 

x 8 ' (u ) (hi du dx 


(29) 


Note that we have used 

-^ = - sin8 (30) 

du 1 c 


and g( 2 ) must be restricted to/= 0 because we are deal- 
ing with 8' (u ) where the curvature term of 8 \u ) 
has already been removed (see Eq. (18)) and added to 
(h . We must mention here that, in new variables 
q 2 (y,T) = q 2 [y(u , u , 0,x),x] . 

The condition g = 0 in Eq. (29) implies that 
u = u (u , u , x) . We will use this result in the inte- 
gration of 8' (m 1 ) in Eq. (29). We get 


4jc(|) 2 (x, x) 


= " C J 


du 


q^HLfy 

r sin 8 


du 2 dx 


g = 0 J u =0 


We have the following results 
l| 


du_ 

du 


cote 


g = 0 


(31) 


(32) 


Osin 8 


du 


cote 


g = 0 


(33) 


f =- r H ^ cot0 


(34) 


= - cote v • t x 8(7) 

du (35) 


where t x is the unit vector along the projection of r on 
the local tangent plane and r J . (sum on i) is the 
Christoffel symbol of second kind in the new coordinate 
system. Note that we have a locally orthogonal frame 
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which gives g( 2 ) = 1 at the origin but because we have a 
curved surface, we get Eq. (34). Therefore, in the rest of 
the algebraic manipulations, we will set g( 2 ) = 1- 

We can now write Eq. (3 1 ) in the following form 
after taking the derivative with respect to m 3 of the inte- 
grand: 

r F cos0 _ „ , cos0 1,2. 

= c — =- t x ■ v 2 q 2 + ~ — du ch 

'-r sin 0 r sin 0 

+ c J c °^ T j" q 1 dir dx 
rsin 0 

v • COS0 

+ c | r q 2 dx 

rsin 0|cos\|/| (36) 

where \|/ is the angle between/ and the edge of the 
open surface described by / = / = 0 . 

In our previous work, we have used dT for du 2 
where dT is the element of the curve of intersection of 
/= 0 with the collapsing sphere g = 0. We have shown 
also that 6,10 


cdxdT _ dl 
sin0 A 

(37) 

dx _ dL 
cos\|/ A 0 

(38) 

where L is the edge of the I -surface 
F = F = 0 and 

described by 

A 0 = |VFxVF| = A A sin 0' 

(39) 

cos0' = N-N' 

(40) 

Zt _ VF _ v-M v r 

\vf\ A 

(41) 

A 2 = 1 +M 2 -2M v cos0 

(42) 

cos0 = v • r = v • fj 

(43) 


Therefore, Eq. (36) can be written as follows: 



We have separated the integrals over the X - surface in 
this equation to simplify the analysis of the singularities. 


Finally, the solution of Eq. (25) was also given by 
Farassat 2 ’ 5 6 as follows: 



We have thus given the full solution of Eq. (19) which 
we refer to as formulation 4. Only Eq. (44) in our analy- 
sis is new. In the next section, we will discuss the 
important question of the singularity of the solution of 
Eq. (19). 

4. A Study of the Singularities of the Solution for 
Supersonic Surfaces 

We will now address the question of the singulari- 
ties in the solution of the FW-H and K equations. There 
has been a general belief among the researchers, the 
authors of this paper included, that the solution of these 
equations lead to nonintegrable singularities for some 
observer space-time variables (x, t). We will show here 
that for a smooth surface, all the singularities of the new 
solution are integrable. We assume that /= 0 is not an 
open surface and thus we only consider the integrals 
over the Z - surface : F(y ; x, f) =fiy, t - r/c) = 0. See 
articles by Farassat, De Bernardis and Myers for the 
analysis of the singularities of the line integrals 11,12 . 
As will be seen below, the analysis of the singularities 
of the surface integrals is very difficult. 
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There are two kinds of singularities in the surface 
integrals: i) sin0 = 0, i. e., 0 = 0° or 180° but 
A 0 , and ii) A = 0 . We will show below that A = 0 
also implies 0 = 0° or 180° . Condition i) means that at 
some source time, the collapsing sphere g = 0 is tangent 
to the surface/= 0 at a point where M n ^ ±1 . Condition 
ii) means that at some source time g = 0 is tangent to 
/= 0 at a point where M n = ±1 . We will prove these 
assertions for appearance of singularity A = 0 below. 
We assume that/= 0 is convex with no saddle points. 

Consider a rotating surface part of which moves at 
supersonic speed. We can write A as follows: 


a constant. We next use the relation 10 

cTL _ c clT dx , 

X sintT ■ • 

where T is the curve of intersection of g = 0 with/ = 0. 
Near the point A, above, dY = Me p, sin0 = b/r , 
where b is the radius of r which is a small circle, and cp 
is the azimuthal angle around T. Thus, we have 

crdtydx (51) 

which means if q(y, x) is continuous, then 


A 2 = (1-M„cos0) 2 +M^sin 2 0. (46) 

Therefore, A = 0 if sin0 = 0 and 1 -M n cos0 = 0 
simultaneously. This means that we must have M n = 1 
and 0 = 0° or M n = -1 and 0 = 180°. The 
geometrical interpretation of these conditions is obvi- 
ous. Since /= 0 is assumed to be moving supersonically 
on part of its surface, there is a curve 'P on/= 0 on 
which M n = 1 or M n = -1 . If at any source time x 0 , 
the collapsing sphere g = x-t + r/ c = 0 is tangent to 
/= 0 at a point on ¥ , then the integrands of the surface 
integrals of formulation 4 are singular for the time 
t 0 = x 0 + r/c . Note that, in general, the curve ¥ is 
time dependent but it is not so for a hovering rotor oper- 
ating at supersonic tip speed. This is the case we study 
here. 

We must study the two conditions of appearance of 
singularities separately. The reason becomes apparent 
below. Let us first write the kinds of integrals we have: 



F> 0 (47) 


r i-23t 

7 1 = C l dZ \ 0 ^2\et d ^ (52) 

X 

is integrable. Therefore, we have shown the integral 


h - J TR^reP 

F = 0 

F> o (53) 

is integrable. 

For I 2 , we use the fact that sin0 = b/r and 
b = Cj jxj to write the integral in the form 


h 


t cadx 

\cM 


2k . 

I r li [ Q2\ et d<P 


(54) 


where a = r(x = 0) . This integral is convergent. 

The study of convergence of 7 3 is very interesting. 
It appears that this integral is not convergent. We 
manipulate the integrand as follows near the condition 
sin0 = 0 : 


E = 


COS0 

rsin0 


h ' V 2#2 


T _ f COS0 „ i , i jt- 

2 1 rAsin0 1?' ret 

F = 0 

F> 0 (48) 


7 3 = J 

F = 0 
F> 0 


' COS0 
3-sin0 1 


• V 2^' 


COS 0 

2 • 2 n qi \ 

r sin 0 J ret 


dL 

X 


(49) 


We will show that all the singularities are integrable. 


Condition i): sin0 = 0. A t. 0 

Near this point, the intersection of g = 0 and/= 0 is 
a circle of radius b. We can show that if this condition 
appears at x = 0, then as x — > 0 , b = Cj fxj where C is 


cos0 a , V 7 

+ — q 2 - —[b t t ■ V 2 q 2 + q 2 ] 

r sin 0 rb 


(55) 


We note that / ■ V 2 ^ 2 = ~dq 2 /db and 


b fj • y 2 q 2 + q 2 = b 


-.2 

2 d _d2 
db 2 


(56) 


and, thus, near the point of tangency of /= 0 and g = 0: 
^2 

( 57 ) 


„ a & q 2 
E ~ - 


r db 2 
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The convergence of / 3 , therefore, depends on the value 
of the following integral when sin0 — > 0 : 


/ 3 = c flj dx J 


2jc 


0 


d q 2 
db 2 . 


r/cp 


x = 0 


(58) 


We can easily show that 
-27t f 


r ^V 


0 




r/cp = [V 2 # 2 ] 


x = 0 


J T = 0 


(59) 


The convergence of / 3 is, thus, guaranteed. We conclude 
that when sin0 -» 0 and A ^ 0 , formulation 4 has only 
removable singularities. 


Condition ii): A = 0 (implies sin0 = 0) 

As we have shown above, when g = 0 is tangent to 
/= 0 at a point where M n = ±1 , we have A = 0 and 
from the tangency condition sin0 = 0 . The first thing 
we study here is the structure of the E - surface near a 
point where A = 0 . We then study the problem of the 
singularities of formulation 4. 

We consider the condition 0 = 0° and M n = 1 
for a hovering rotor. Figure 1 shows the tangent plane T 
to a point on V looking edgewise at the moment of tan- 
gency of the collapsing sphere with/= 0 at the point A. 
Note that M n = n = r at this moment so that the 
observer is in the plane shown at the center of the col- 
lapsing sphere g = 0. We assume x 0 = 0 and r = a at 
the moment of tangency. Since M n = 1 at A, we have 
ship. = 1/M where p is the angle that T makes with 
Mach number vector M shown in Fig. 1 . We now con- 
sider the plane T in motion for |x| < £ where £ > 0 is a 
small number. In the frame fixed to T with origin at A as 
shown in Fig. 2, the curve of intersection of g = 0 with 
T, which is a circle is given by the relation 


(x-pv a x) 2 + (y + yV R x) 2 = v 2 j\ 2 (60) 


Where V a = aco, V R = R (0 and P and y are defined 
in fig. 2. We use R for the distance from A to the axis of 
rotation. The center of the circle is at the point 
(x ,y) = ($V a x, -yV R x) and its radius is C a P[x| . 
It is clear that as |x| — > 0 , the E - surface looks like a 
vertex of a cone and has no tangent at the vertex coin- 
ciding with A. The condition of A = 0 is thus 
equivalent to the E- surface becoming pointed and 
having no tangent plane at the point A. 

Let us see what the intersection of g = 0 and T will 
appear to an observer on the tangent plane T. Figure 3 


shows the envelope of the circles of Eq. (60). We reject 
the part of the envelope in the region M < 1 because it 
would imply multiple emission from subsonic region 
which is impossible. This means that we have no inter- 
section of g = 0 and T for x < 0. This figure clearly 
shows that the E - surface looks like a cone near A. 



Plane 

Fig. 1 . The geometric condition for the appearance of 
A = 0: 0 = 0°, M n = 1 , i.e., M n = r = n . 



Fig. 2. The coordinate system used to study the intersec- 
tion of g = 0 with/= 0 near the condition A = 0. 


For convergence of I\, the analysis of condition i) 
applies exactly so that 7) has removable singularity. The 
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y tan a = (R tan ja)/a 



z-axis out of the plane 

Fig. 3. The shaded area is the trace of the intersection of 
g = 0 in the xy-plane of figure 2 near the condition 
A = 0. It is assumed that x = 0 for A = 0. 

study of the convergence of I 2 is different since from 
Eq. (60) we see that as A— >0, we have b = Cjt 
where C\ is a constant. However, we note that for both 
the FW-H and K equations, q 2 is proportional to 

2 

M n — 1 . We can easily show that as A — > 0 , 

2 

M n - 1 = C 2 x where C 2 is another constant. There- 
fore, the convergence of I 2 is guaranteed because near 
A = 0 , we can write I 2 as 



2 

where q' 2 = q 2 /(M n - 1 ) . As a matter of fact, / 2 , is 
better behaved for condition ii) than condition i)! 

For/ 3 , the convergence study of condition i) applies 
exactly for condition ii). Therefore formulation 4 has 
only removable singularities for condition ii) also. 

We conclude that for a smooth surface, the solution 
of FW-H and the K equations as given here have only 
integrable (removable) singularities for a supersonic 
surface f = 0. The solution of the K equation is, of 
course, known as the Kirchhoff formula for supersonic 
surfaces. 


Remark 1 . The solutions of the FW-H and the K 
equation here are valid for all range of the surface speed. 
But we do not recommend to use the present results for 
subsonic surfaces since much more efficient solutions 
for numerical method are available 2,5,10 . 

Remark 2 . It can be shown that had we not added 
the surface terms from the quadrupole source term of 
the FW-H equation to the thickness and loading source 
terms, the resulting solution would be singular when the 
condition A = 0 appears. The acoustic pressure signa- 
ture will have a logarithmic singularity which will 
appear as an infinite pulse. Our analysis shows that 
when all the surface sources from thickness, loading and 
quadrupole terms are included in the analysis, there is 
no infinite singularities in the acoustic field. 

5. Concluding Remarks 

We have given the closed form solution of the 
FW-H and the K equations for an open surface. 
Although these solutions are valid for all range of Mach 
numbers, we recommend them for the supersonic 
motion of the surface because of their complexity. We 
have shown that for a smooth surface, the singularities 
of the solutions of both of these equations are integrable. 
The nature of these singularities is explained in this 
paper. It is very interesting to note that for the FW-H 
equation, the thickness and loading source terms alone 
have nonintegrable singularities in the solution. How- 
ever, the addition of the surface source terms from the 
quadrupole source term removes this singularity. This 
is, of course, expected on intuitive grounds. 

We hope that the present work gives further impe- 
tus to numerical applications of our results in high speed 
rotating blade noise prediction. The closed form analytic 
results of this paper open up two other areas of applica- 
tion which could help aeroacousticians in their endeavor 
to reduce the noise of aeronautical machines. These 
areas are: i) qualitative analysis of noise generation 
mechanisms by the analytic study of the appropriate 
integrals in our solutions of the FW-H and K equations, 
and ii ) approximate analysis of the radiation field from 
ducted fan inlet and exhaust and other openings that 
radiate sound to an infinite medium. The analysis is sim- 
ilar to the use of the conventional Kirchhoff formula for 
the study of diffraction by an aperture. 

Much work is ahead of us in the order of magnitude 
analysis of the terms in the solution of the FW-H and the 
K equation. Can a deformable body /= 0 be used to 
control noise radiation at high speed? Can we design a 
rigid body with desirable noise radiation property in a 
given direction? 
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Errata For AIAA 98-2375: A Study of 
Supersonic Sources- The Ffowcs Williams- 
Hawkings Equation and the Kirchhoff 

Formula 


F. Farassat, Kenneth S. Brentner and M. H. Dunn 

There is an error in Eq. (30) of this paper presented at the 4th AIAA/CEAS 
Aeroacoustics Conference in Toulouse, France resulting in a number of 
changes in the solution of the wave Eq.(28). Equation (30) must be corrected 
to 


dg_ 

di / 1 


1 1 3 

-[sin0 + K,(M cosO -u sinO)] 
c 1 


(30) 


where Kj is the normal curvature of the surface / = 0 in the direction of t j 
which is the unit vector in the direction of the projection of the radiation vec- 
tor r on the local tangent plane. This relation is valid to the first order in u 

3 

and u . Equation (29) must be written as follows: 


4jc<|> 2 (x, t) = J 


q 2 (y^) 


3 


H(f) 8' (w )8 (g) g, r..du du du dx 

r A a/~(2) 


1,2,3 


h4 8 <2) H{f) 


dg/du 


(29) 


3 2 3 

5 ' (u ) dll du dX 


Jg = 0 


Equation (31) becomes 


4tz<\> 2 {x, t) = - 


(k 


q 2 H (f)Js^) 

r\dg/du l \ 


du 2 dx 


(31) 


J g = 0 J u =0 


Errata For AIAA 98-2375: A Study of Supersonic Sources- The Ffowcs Williams-Hawkings Equation and the Kirchhoff Formula 



Initially we assume that the collapsing sphere is not tangent to the panel 
/ = 0, / > 0 as it crosses the panel. It can be shown that 


" a 

- 

1 

1 ^ 

COS0 K 1 

law 3 

L r 

dg/du 

s = °V = o 

_ r 2 (sin0) 3 r(sin0) 3 _ 


Then, using the above result and after some algebraic manipulations, Eq. (36) 
becomes 


4 7 T 4 > ^ ( -XT, t) 


C J 1 

j r 


COS 0 . V 7 K 1 

,Ti:'r v 2«2 -sm«2 

Lsin 0 


2 

dll dX 


V • ty COt 0 
r|cos\|/| ^2 


dX 


(36) 


Equation (44) becomes 


1 fj ^ 

= J -U'oier, V 2 </ 2 - K,« 2 l ref x 




V • ty COt 0 

a7 


q 2 


dL 


J ret 


(44) 


If, however, the collapsing sphere leaves the panel tangentially, another line 
integral similar to that in eq.(44) around the edge of a hole enclosing the 
point of tangency must be added. The limit of this line integral, as the maxi- 
mum diameter of the hole goes to zero, adds the following term to Eq.(44) 
which is not in the paper: 


/ = 


4^2 f 1 ^ 2 sig[k((p)] 


if; 


1 -M r J 0 k r -k( cp) 


chp 


Jr 


Here we have defined k f = l/r, sig(.) is the signum function and k(ty) is 

the normal curvature of the surface / = 0 at the point of tangency T as a 
function of the azimuthal angle cp . The sign convention for the curvature is 
based on assuming n T = r T , i.e., /c(cp) > 0 if the center of curvature is on 


Errata For AIAA 98-2375: A Study of Supersonic Sources- The Ffowcs Williams-Hawkings Equation and the Kirchhoff Formula 



the same side that n T points into. We assume that the surface / = 0 has 

nonnegative Gaussian curvature everywhere. This is not a severe restriction 
on the surface / = 0 since, in general, one avoids a surface with saddle 
points. Note that at the point of tangency T we have M r = ±A/ ; . Since t/ 0 
2 

has a factor of M n - 1 , the above equation is not singular even when the col- 
lapsing equation leaves the panel at a point where there is a Doppler singular- 

2 2 

ity. Using Euler’s formula &(cp) = K/coscp) + K 2 (sin(p) , where K, and 

k ~2 are the principal curvatures at the point of tangency, the above integral 

can be integrated in closed form with respect to (p . Under some conditions, 
e.g., when 0 < k 2 < k < Kq , this integral must be interpreted as the principal 

value integral. The limit for a flat point or a cylindrical point can also be 
obtained. The above integral appears in geometrical acoustics and geometric 
diffraction theory. The full discussion of this point as well as the verification 

of the final results will be published soon 1,2 . The conclusions of the paper are 
correct. Specifically, we claim that we have presented the simplest possible 
formula (designated Formulation 4) for prediction of the noise from high 
speed (transonic and supersonic) moving surfaces . The discussion of singu- 
larities in the paper must be changed in light of the above results. Elowever, 
the corrections given here improve the behavior of the integrals in Formula- 
tion 4 at the singularities which are all removable. 

We take this opportunity to bring to the attention of the readers the following 
misprints in this paper: 

1 . Equation (29), p 4, replace/ with f 

2. Page 4, second column 6'//) and 5' (z/) must be replaced everywhere 

3 3 

by 6'(zz ) and 5' (z/ ), respectively 

3. Equation (31), p 4, replace by 

4. The left sides of Eqs. (34) and (35) are evaluated at g = 0 
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5. Equation (38), p 5, replace dx with cdx 


6. Equation (40), p 5, replace N' with N 

7. Equation (43), p 5, replace t) ■ ty with \) ■ sinO 

8. First paragraph of Sec. 4, p 5, third line, ‘brief’ must be replaced by 
‘belief’ 

The authors thank Professor Mark Farris of Midwestern State University, 
Wichita Falls, Texas who pointed out the error in Eq.(30) and independently 
verified our results in detail. He was the 1998 ASEE Summer Faculty Fellow 
at NASA Langley Research Center working with F. Farassat. 

References 

1. F. Farassat, Kenneth S. Brentner and M. H. Dunn: A New Formula for the 
Study of the Noise From High Speed Surfaces- The Ffowcs Williams-Hawk- 
ings Equation and the Kirchhoff Formula, submitted to the Journal of Sound 
and Vibration 

2. F. Farassat and Mark Farris: Verification and Analysis of Formulation 4 of 
Langley for the Study of Noise Radiation From High Speed Surfaces, sub- 
mitted to the 5th AIAA/CEAS Aeroacoustics Conference, 10-12 May, 1999 


Errata For AIAA 98-2375: A Study of Supersonic Sources- The Ffowcs Williams-Hawkings Equation and the Kirchhoff Formula 



